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Abstract 



We address the problem of determining whether or not a harmonic oscillator has 
been perturbed by an external force. Quantum detection and estimation theory has 
been used in devising optimum measurement schemes. Detection probability has 
been evaluated for different initial state preparations of oscillator. The correspond- 
ing lower bounds on minimum detectable perturbation intensity has been evaluated 
and a general bound for random phase perturbation has been also induced. 



53 ■ 1 Introduction 

a* 

The harmonic oscillator is a relatively simple model, which is widely utilized in 
many fields of physics. Indeed, it provides a satisfactory description of a large 
number of very different physical systems. This is true also in a quantum mechanical 
framework, where the harmonic oscillator plays a crucial role. Its spectrum of 
eigenvalues, in fact, is infinite, discrete and bounded from below, thus representing 
a paradigm for any bounded oscillating system. 

Interesting physical features often comes with perturbation to harmonic be- 
haviour, which are to be revealed from measurement performed on the system. 
This is the case, as an example, of oscillating electronic circuits [[jj or of some large 
mass viewed as a gravitational antenna [^, [| |J. Also, a single mode radiation field 
is modeled on harmonic oscillator, and many optical devices act as driving terms in 
the dynamical evolution [||. 

Quantum mechanically, different measurements lead to different informations 
about a physical system, as each observable shows up only an aspect of the state 
under examination ffl, ffl, H. Therefore, it is a matter of interest to analyze the 
different measurement processes, in order to find an optimized measurement scheme, 
which is capable to reveal perturbations as weak as possible. This project is a matter 
of quantum detection and estimation theory which regards the very general 

problem of extracting information on a physical system from measurement or a set 
of measurements. 
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In the present paper we address the binary decision theory for driven quantum 
harmonic oscillator. Let us consider an oscillator which is free to follow harmonic 
evolution, and possibly subjected to an external driving signal. After a fixed time 
r we perform a measurement on the system, in order to check whether or not the 
oscillator has been perturbed. We are going to deal with two questions: first, which 
is the best measurement one can perform, in order to reveal perturbations as weak as 
possible with the minimum probability of error ? And second, which is the minimum 
detectable perturbation intensity, depending on the initial state preparation of the 
oscillator ? 

We do not concern to any specific measurement device and we do not discuss 
the feasibility of optimized measurement. Rather, we attempt to derive a ultimate 
quantum limit on the detectable intensity of a perturbation, which depends only on 
the initial quantum state of the oscillator. 

The paper will be organized as follows. In Section |2| we set the proper quan- 
tum measurement theory framework and illustrate the Neyman-Pearson strategy for 
binary decision. In Section || we consider different initial preparation states for the 
harmonic oscillator and derive the corresponding minimum detectable perturbation 
intensity. Section ^| closes the paper with some concluding remarks. 



2 Quantum Detection Theory 
2.1 Quantum measurements 

In a quantum mechanical framework any measurement apparatus is a device, at 
least a mathematical one, which turns each quantum state p into a probability 
density distribution 

dP[p]: P — >dP[p](x) xeX, (1) 

where X is some measurable space, where the possible outcomes of the measurement 
lie. This can be the real Borel set or a subspace of it. The measurement map is 
provided by trace operation 

dP[p](x)=Tr{pdp(x)} , (2) 
which assures propagation of convex linear combinations from density operators to- 



ward probabilities. The above formula is the Born's statistical rule [12|. It contains 
the whole probabilistic structure of quantum mechanics [jl3| . The Born rule leads to 
a genuine probability density distribution if the operator d(i{x) satisfies the axioms 



for a probability operator measure (POM) namely it is nonnegative 

dpl(x) > , (3) 
and it provides a resolution of identity on the set of possible outcomes 

dp{x) = i . (4) 

x 

Eq. (||) guarantees the probability density in Eq. (|2j) to be normalized, whereas 
positiveness of dp(x) also assures it is selfadjoint. 
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Spectral, orthogonal resolution dE(z) of a selfadjoint operator 

Z = J z dE{z) dE{z) = \z){z\ dz 

Z\z)=z\z) (z\z'} = 5 z (z - z') , (5) 

provides a projection valued measure (PVM) which belongs to the class of POM. 
However, this is not the most general example. Also nonorthogonal projectors or 
overcomplete sets provides POM, namely available measurement schem e [jl5] , [l7[|. 

It is worth noting that a POM is necessarily (Naimark Theorem) [lg] a partial 
trace of a PVM coming from a selfadjoint operator defined on a larger Hilbert space. 
The latter can be thought as the whole Hilbert space describing both, the system 
under examinations and the measurement apparatus []l9| . As we are not going to 
deal with physical implementation of measurement we can restrict our attention on 
the Hilbert space of the examined system only. Thus, any measurement is properly 
described by a POM. 



2.2 Driven Harmonic Oscillator 

The quantum mechanical description of harmonic oscillator is based on annihilation 
a and creation operators 



mm / 1 

which form the number operator N = a^a. Due to commutation relation [a, a^] = 
1, multiple applications of to the vacuum state leads to the Fock basis |n) = 
(n\)-^ 2 a^ n \0) which span the whole Hilbert space representing the possible levels 
of excitation for the oscillator. Number states represents also the eigenstates of the 
number operator, whose spectrum coincides with the set of the natural numbers 
J\f = 0,1,.... The eigenstates of annihilation operator a\a) = a\a) constitute the 
overcomplete set of coherent states, a being the complex amplitude of the harmonic 
oscillations. Coherent states can also be obtained from the vacuum by the action 
of displacement operator D(a) = exp{<W — aa} (2(| 



\a) = D(a) 

b(a)b{a') = l)(a + a / )exp{Im(aa / )} . (7) 
Let us consider the Hamiltonian (natural unit % = 1) 

H = -t^ + ^ v+Fit)x - (8) 

It describes a classical harmonic oscillator of mass m and frequency uj subjected to 
a time dependent driving force. In terms of annihilation and creation operator the 



Hamiltonian could be written as |21] 



Fit) 

H = ojo)a H — j=J=(a) + a) . (9) 

\/2muj 
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Finally, we adopt interaction (Dirac) picture to obtain 

Hj = ^22=(at + a). (10) 

Let now consider the initial state of the oscillator to be pq. If no driving force is 
present the final state, after a fixed evolution time r, is still po- Otherwise, we have 

p 1 = Up tf , (11) 

where the evolution operator U = exp {iH/r} could be written as a displacement 

U = D{z) = exp {za) - za} , (12) 



operator [22] 



where 



7 = / dte iujt F(t) . (13) 



y/2mu> ' Jo 

The quantity z represents the complex amplitude of the driving signal, whereas \z\ 2 
denotes the energy intensity of the perturbation, expressed in unit of the oscillator 
quanta u. 

2.3 Neyman-Pearson Strategy for Binary Decision 

Our goal is to determine whether or not the system has been perturbed. For this 
purpose we adopt a detection scheme as in Fig. |l|. After the initial preparation the 
harmonic oscillator is left free to evolve for a fixed time r. Then, some kind of mea- 
surement dp(x) is performed. Starting from the outcomes of such a measurement 
we have to infer which is the state of the system, in order to discriminate between 
the following two hypothesis: 

TCo : No perturbation has been occurred during the time interval r, true if we infer 

H\\ The system has been perturbed during the time interval t, true if we infer p\. 

We denote by Poi the probability of wrong inference, namely that one of inferring 
Ji\ when 7io is true. In hypothesis testing formulation this is usually referred to as 



false alarm probability |23|]. Conversely, we denote by P\\ the detection probability, 
that is the probability of inferring 7i\ when it is actually true. 

Now, which is the best measurement to discriminate between po and p\ ? 

If these two states are mutually orthogonal the problem has a trivial solution. 
It is a matter of measuring the observable for which po and p\ are eigenstates. 
However this is not our case, as displacing a state of the harmonic oscillator leads 
to a different kind of state. Only coherent states maintain their characteristic under 
displacement 

t){z)\a) = D(z)D(a)\0) = exp(Im[az])\a + z) . (14) 

However, coherent states constitute a nonorthogonal, overcomplete set by them- 
selves. Thus, the above procedure cannot be applied in the present case, even for 
special initial states of the oscillator. 
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figure 1: outline of the detection scheme. the system is initially prepared 
in some fixed state and then it is left free to evolve. during the free 
evolution it could be subjected to an external driving signal. after a 
fixed time interval t we perform some kind of measuremen on the system. 
From the outcome of such a measurement we have to infer which is the 
state of the system, in order to discriminate between the perturbation 
hypothesis and the null hypothesis. 



In the following we consider nonorthogonal po and p\ and we focus our attention 
on oscillator initially prepared in a pure state po = \ipo)(ipo\. As it can easily checked 
from Eq. flIT| ) this means that also the perturbed state is a pure state p\ = 

The optimization problem can be analytically solved, for pure states, by adopt- 
ing, the Neyman- Pearson criteria for binary decision p4}| . The latter reads as follows. 
First, we have to fix a value for the false alarm probability Pq\. Then, we have to 
find the measurement strategy dfi{x) which maximizes the detection probability P\\. 
As a general definition, each measurement strategy which maximizes the detection 
probability P\\ for a fixed value of false alarm probability Poi is considered as a 
Neyman-Pearson optimized detection for binary hypothesis testing. It was shown 



by Helstrom ||] and Holevo |1C] that this very general problem could be reduced to 



solving the eigenvalue problem for the operator 

dp(x\X) = pi - Xp , (15) 

which represents the optimized measurement scheme. In general it is a POM rather 
than a PVM. Nevertheless when, as it is here the case, the two signals are linearly 



independent it has been proved by Kennedy [25, [2(| that the optimum detection is 
described by a PVM. The parameter A is a Lagrange multiplier. Different values of 
A correspond to different values of the false alarm probability, namely to a different 
Neyman-Pearson strategies. 

Once the eigenvalues problem for dfi(x\X) has been solved it results that only 
positive eigenvectors contribute to the detection probability P\\ |], |27|, p8[ . Thus the 
decision strategy is transparent: after a measurement of the quantity dp{x\\) if the 
outcome is positive we infer perturbation hypothesis H\ is true. Conversely, we infer 
null hypothesis TCq when obtaining negative outcome. By expanding the eigenstates 
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of dfi(x\X) in terms of \tjjo) an d Lagrange multiplier A can be eliminated from 
the expression of detection probability which results 



P 
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V7^+V(l-P i)(l 



< Pqi < K 



« < Poi < 1 



(16) 



In Eq. ( |16| ) k denotes the square modulus of the overlap between perturbed and 
unperturbed state of the harmonic oscillator, in formula k = \0[tPq, z]\ 2 , where 



0[ih,z] = (ih\ih) = (^\D(z)\iM 



(17) 



The overlap depends both on the initial state and on the perturbation amplitude. 
In the next Section we evaluate the quantity in Eq. ( R% for relevant kinds of initial 
state. 
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Figure 2: THREE DIMENSIONAL PLOT OF THE DETECTION PROBABILITY P u FOR 
NEYMAN-PEARSON OPTIMIZED DETECTIONS AS A FUNCTION OF THE FALSE ALARM 
PROBABILITY P 01 AND THE OVERLAP K BETWEEN THE INITIAL STATE AND THE PER- 
TURBED ONE. 



In Fig. U we report the detection probability of optimized detection strategies 
as a function of the false alarm probability Pqi and the overlap strength parameter 
k. It is obvious that if the overlap is small, it is easy to discriminate between the 
two states. Thus, it is possible to obtain strategies with large detection probability 
without paying the price of an also large false alarm probability. On the contrary, if 
the overlap becomes appreciable it is difficult to discriminate the states. In the limit 
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of complete overlap the perturbed and the unperturbed states become indistinguish- 
able. Detection probability is now equal to false alarm probability and the decision 
strategy is just a matter of guessing after each random measurement outcome. 

Choosing a value for the false alarm probability is a matter of convenience, de- 
pending on the specific problem this approach would be applied. The maximum 
tolerable value for Poi increases with the expected number of measurement out- 
comes, and conversely a very low rate detection scheme needs a very small false 
alarm probability. 



3 Quantum Binary Decision for Harmonic Os- 
cillator 

Once an acceptable value of false alarm probability has been fixed and the oscillator 
has been prepared in some initial state \ipo), the detection probability P\\ depends 
only on the perturbation intensity. A wise inference could be performed only when 
Pn( z ) > 1/2, as only in this case the record X = {xo,x±, ...,xn} of experimental 
data contain usable information. Thus, the threshold value P\i{z) = 1/2 defines 
the minimum detectable perturbation for the so-prepared oscillator plus detector 
system. We will consider the intensity of the minimum detectable perturbation 
as the relevant parameter and we denote it by \z m i n \ 2 = M.. In the following 
Subsections we study the behaviour of k and M for different initial preparation 



3.1 Coherent States 

For the oscillator prepared in a coherent state \a) the overlap is given by 

0[a, z] = (a\D(z)\a) = exp{-^|z| 2 } expjza - az} , (18) 

and thus the overlap strength does not depend on the amplitude of the prepared 
coherent state 

k = exp{-|z| 2 } . (19) 

For zero false alarm probability the detection probability is given by P\i{z) = 
1 — exp{ — |z| 2 }. When a small false alarm probability is set (Poi < K ) the minimum 
detectable perturbation intensity is obtained by the inversion of the formula 



V^^+^/(l-Poi)(l-«) 



(20) 



that is, 



^MwW (21) 



The minimum detectable intensity is independent on the initial coherent amplitude. 
Thus coherent states provide a stable oscillating system, however also difficult to 
monitor in its fluctuations. 
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3.2 Squeezed States 



Uncertainty principle set a lower bound for the product of fluctuations for two 
conjugated quantity. This implies a degree of freedom, namely that one can ar- 
bitrarily reduce the fluctuations in some variable upon increasing the fluctuations 
in the conjugated one. Indeed, squeezed states of the harmonic oscillator have 
been introduced as minimum uncertainty state for amplitude quadrature operators 
Xtp oc a'e ilp + ae~ lip with phase dependent fluctuations. 

Squeezed state can be obtained by the coherent displacement of squeezed vacuum 
\a, Q = D(a)\O,0 |29| ) [30[| . The latter is obtained from the vacuum by the action 
of squeezing operator |0, () = S(()\0), where 



S(C) = exp^ Ca t2 -Ca 



(22) 



and C = rexp{i2ip}, with r real. Squeezing a state implies the introduction of 
some energy. The mean excitation number of a squeezed vacuum is given by n sq = 
(C,0|iV|0,C) = sinh 2 r. 

As we have seen just above coherent amplitude does not cause any effect when 
subjected to displacement action. Therefore, we restrict our attention to squeezed 
vacuum |0, £) which shows all the interesting phase dependent features related to 
squeezing. We also consider, for simplicity, a squeezed vacuum with squeezing phase 
equal to zero tp = 0. 

The overlap is given by 



C[r, z] = exp | — -|z| 2 cosh 2r — sinh 2r cos 2 <p j 



(23) 



where <p = arg(z) is the phase of the perturbation. In Fig. || we report the overlap 
(|23| ) for a unit intensity perturbation \z\ 2 = 1 as a function of the squeezing param- 
eter r and the perturbation phase (p. In the two limiting cases ip = 0,7r/2 we have, 
for the overlap strength 



kq = exp 



-2\z\ 



+ 



i n/2 = exp • 



1 

2 ~ 
2 [2n 



(n sq + l)n sq 



sq 



(24) 



(25) 



However, the perturbation phase is reasonably random, or unknown. Thus, a rele- 
vant parameter to be considered is also the phase averaged overlap strength which 
is defined by 

r n dp 



2tt 



0[r, Me^] 



(26) 



Inserting Eq. (p3|) in Eq. (pBf) leads to 



exp 



2n sq + 1 - %J (n sq + 1) 

2 



(27) 



with Iq(x) denoting a modified Bessel function of the first kind |31]. 



S 
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Figure 3: CONTOUR PLOT OF THE OVERLAP FOR THE OSCILLATOR INITIALLY PRE- 
PARED IN A SQUEEZED VACUUM. A PERTURBATION WITH UNIT INTENSITY IS CONSID- 
ERED AND THE OVERLAP IS REPORTED AS A FUNCTION OF THE SQUEEZING PARAME- 
TER r AND THE PHASE (f OF THE PERTURBATION. 



The minimum perturbation equation (|20|) can be analytically solved for a fixed 
value of perturbation phase. We obtain 



M = log 



1 + ^01(1-^01) 



exp{2r} 



(28) 



where e 



2r 



2n sq + l + 2^J {n sq + l)n sq , and 

2 



M 



n/2 



log 



1 + ^01(1-^ 



01 J 



2n sq + 1 



(29) 



From Eqs. ( |28| ) and ( p9| ) is apparent the strong effect of phase matching. When 
both phases, the perturbation one and the squeezing one, have the same value the 
overlap is strongly enhanced and thus the minimum detectable intensity increase 
(roughly linearly) with the increasing of the squeezing energy. On the contrary when 
the two phases are maximally mismatched, the overlap decreases with increasing 



9 



energy of the initial states. Thus, the system becomes more and more sensitive to 
perturbation and minimum detectable intensity shows an inverse scaling with the 
initial preparation energy. 




Figure 4: THE MINIMUM DETECTABLE PERTURBATION INTENSITY M FOR THE OSCIL- 
LATOR INITIALLY PREPARED IN A SQUEEZED VACUUM AND RANDOM PHASE PERTURBA- 
TION. The behaviour of M is reported as a function of the mean excitation 

NUMBER n sq = sinhV OF THE SQUEEZED VACUUM FOR DIFFERENT VALUES OF THE 
FALSE ALARM PROBABILITY P 1- THE CURVES ARE CLEARLY DISTINGUISHABLE IN 
THE REGION n sg ~ 1 WHERE WE HAVE, FROM TOP TO BOTTOM, THE BEHAVIOUR FOR 

P i = 0.00, 0.01, 0.02, 0.05 respectively. 

In the case of random (unknown) phase we have not been able to solve analyti- 
cally the perturbation equation. We solved it numerically. In Fig. |I| we report the 
behaviour of M as a function of n sq for various values of the false alarm probability. 
For 

n sq ^ 10 the numerical results are very well interpolated by the formula 

, d n sq >10 ( 2 \ 1 . . 

M ~ lOg ; , . 30 

\l + ^P m {l-P m ) n sq 
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3.3 Number States 

The overlap of a number state with its displaced version is a real number, thus the 
overlap strength k is just the square of the overlap 0[n, z\. We have 

k = e W {-\z\ 2 }L 2 n (\z\ 2 ) , (31) 

where L n (x) denotes a Laguerre polynomials pip . Number states are a phase in- 
sensitive kind of states. Therefore, preparing the harmonic oscillator in such a way 
is equivalent to a phase averaging by default. 
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Figure 5: THREE DIMENSIONAL PLOT OF THE OVERLAP STRENGTH K FOR THE OS- 
CILLATOR INITIALLY PREPARED IN A NUMBER STATE. THE OVERLAP STRENGTH IS 
REPORTED AS A FUNCTION OF THE NUMBER n AND THE PERTURBATION INTENSITY 

\z\ 2 . 

In Fig. H we report the overlap strength as a function of the perturbation 
intensity and the excitation number of harmonic oscillator. 
The minimum perturbation equation (|20| ) reads as follows 

^=exp{-M}Ll(M). (32) 
It could be numerically solved. Minimum detectable intensity scales as 

M ^ - , (33) 
n 
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with the proportionality constant depending on the value of the false alarm proba- 
bility. Roughly we have 

A » 0.3 - \p Q i . (34) 



3.4 Superposition of Coherent States 

We end this section by dealing with superposition states. We consider an analyt- 
ically solvable case which is provided by superpositions of coherent states. Let us 
introduce the two set of states expressed by 



1 



2^1±exp{-2H 2 } 



la) ± 



a» 



(35) 



where |a) denotes a coherent state. These states are known also as even and odd 
Schroedinger cats []32| as they are superposition states containing only even and odd 
number components respectively. The evaluation of the overlap can be carried out 
by means of the operatorial relations 



D(-a)D(z)D(a) 
D(a)D(z)D(a) 



exp{za — az}D(z) 
D(z + 2a) . 



(36) 



We consider for simplicity a as a real number, thus we obtain 

exp{-±|z| 2 } 



0[a,±,z) 



1 ± exp{-2a 2 } 



cos (2a|z| sin p) 



±e 



-2a 1 



cosh (2a|z| cos p) 



(37) 



being ip the perturbation phase. After some calculations we arrive at the overlap 
strength for the fixed value tp = 0, tt/2 and for the phase averaged case. We have 



exp{ — \z\ 2 } 



(1 ± 2exp{-2a 2 }) 2 



1 db exp{— 2a 2 } cosh (2a|z| 



(38) 



V/2 



exp{- 



4a 2 } 



;id=2exp{-2a 2 }) 2 
1 ± exp{2a 2 } cos (2a|z| 



(39) 



exp{ — \z\ 2 } 



(1 ± 2exp{-2a 2 }) 2 



J (2a\z\) ± exp{-2a 2 }/ (2a\z\ 



(40) 



being Jq{x) and Iq(x) the Bessel function and the modified Bessel function of the 
first kind Notice that the mean excitation numbers for the superposition states 
are given by 

|2 l T e- 2 H 2 



n± 



a 



l± e - 2 H s 



(41) 
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The strong effect of phase matching is again apparent. The minimum perturbation 
equation can be solved in the asymptotic region of large excitation numbers, leading 
to the energy scaling 



In the case of random phase perturbation the minimum detectable intensity shows, 
at least asymptotically, an inverse scaling relative to the mean excitation number 
of the initial state. The same behaviour we have observed for squeezed vacuum 
and number state initial preparation, and this seems to indicate a general bound 
for detectability of perturbations. Actually, the minimum detectable intensity M. 
becomes almost independent on the initial preparation in the limit of high excitation 
numbers. 



Quantum detection and estimation theory has been applied in binary hypothe- 
sis testing, regarding possible perturbations on harmonic behaviour of a physical 
system. The action of an external driving signal is described by a displacement 
operator, thus excluding the possibility that the perturbed state of the oscillator 
could be orthogonal to that has been initially prepared. 

The detection probability has been evaluated, for different initial preparation 
of the oscillator, as a function of the perturbation intensity and the initial prepa- 
ration energy. Minimum detectable perturbation intensities have been also evalu- 
ated, which represent the ultimate quantum limit in detecting a perturbation for 
fixed initial preparation of the harmonic oscillator. The lower bounds on detectable 
perturbation show a strong dependence on the phase matching. In-phase perturba- 
tions could be effectively detected only for weakly excited oscillators, as minimum 
detectable intensity linearly increases with initial energy. On the contrary, out-of- 
phase perturbations are easily detected also for high excitations. In the realistic 
case of random phase perturbation, the minimum detectable perturbation intensity 
seems to become independent on the initial preparation, at least in the asymptotic 
regime of large initial energy. This suggests that the inverse scaling M oc n _1 
relative to the initial mean energy could be a general bound. 

An ultimate quantum limit for detection of perturbations could be defined inde- 
pendently from initial preparation, upon a further optimization over all the possible 
quantum states of the oscillator. 
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